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Gate-controlled Kondo screening in graphene: 
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Abstract. - Magnetic impurities in neutral graphene provide a realization of the pseudogap 
Kondo model, which displays a quantum phase transition between phases with screened and 
unscreened impurity moment. Here, we present a detailed study of the pseudogap Kondo model 
with finite chemical potential [i. While carrier doping restores conventional Kondo screening at 
lowest energies, properties of the quantum critical fixed point turn out to influence the behavior 
over a large parameter range. Most importantly, the Kondo temperature Tk shows an extreme 
asymmetry between electron and hole doping. At criticality, depending on the sign of /i, 7k follows 
either the scaling prediction Tk oc with a universal prefactor, or Tk oc \/j,\ x with x « 2.6. This 
asymmetry between electron and hole doping extends well outside the quantum critical regime 
and also implies a qualitative difference in the shape of the tunneling spectra for both signs of \x. 
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Introduction. — Electrons in graphene provide an es- 
sentially perfect realization of two-dimensional (2d) Dirac 
fcrmions [1-3]. Scanning tunneling microscopy (STM) al- 
lows to locally study their electronic structure, includ- 
ing the physics of defects and impurity atoms. Recently, 
Kondo resonances for Co atoms adsorbed on heavily doped 
graphene have been reported [4]. On the theory side, the 
study of magnetic impurities coupled to 2d Dirac electrons 
started in the 1990s in the context of d-wave superconduc- 
tors [5-7]. For a spin- 1/2 impurity coupled to electrons 
with a density of states (DOS) obeying p(u>) oc \u>\ r - the 
so-called pseudogap Kondo problem (with r = 1 for d-wave 
superconductors and graphene) - it has been established 
that Kondo screening exists for low temperatures T only 
if the Kondo coupling J is larger than a critical value J c . 
The properties of the resulting quantum phase transition 
(QPT) have been studied extensively [5-11]. Interestingly, 
the fixed-point structure changes as function of the bath 
exponent r [7]. In particular, for r > r* 0.375 the criti- 
cal fixed point is particle-hole (p-h) asymmetric, and r = 1 
plays the role of an upper critical dimension [10]. 

For an impurity adsorbed on graphene, its position rel- 
ative to the underlying carbon honeycomb lattice and the 
character of its magnetic orbitals determine the electronic 
hybridization paths. As a result, orbital degrees of free- 
dom are crucial for Kondo screening if the impurity is 
located at a high-symmetry position [12]. For the case 



of Co impurities, the most interesting situation has been 
identified as that of a Co atom with spin 1/2 sitting in the 
center of a C hexagon. Here, a spin-orbit-driven crossover 
from SU(4) Kondo physics at high energies to conven- 
tional SU(2) Kondo physics at low energies is predicted. 
Hence, the low-temperature behavior is that of a standard 
Kondo model with a pseudogapped bath DOS. Moreover, 
the Kondo coupling J for this model has been estimated 
to be rather close to J c , corresponding to the QPT be- 
tween screened and unscreened impurity phases in neutral 
graphene [12]. 

This opens the exciting possibility to observe Kondo 
criticality [13] in graphene. Even if a tuning of the cou- 
pling constant exactly to its critical value is impossible, 
quantum critical signatures can still be expected at en- 
ergies above some crossover scale. The most important 
new aspect (as compared to d-wave superconductors) is 
the possibility of doping via a gate voltage, which moves 
the chemical potential fi away from the DOS minimum. 

In this paper, we discuss near-critical Kondo physics in 
the pseudogap Kondo model with finite /i, using a combi- 
nation of analytical and numerical renormalization-group 
(RG) methods, with an eye towards possible experiments. 
As detailed below, we find an extremely strong asymmetry 
in the impurity properties between electron and hole dop- 
ing. This is less related to the p-h asymmetry of the hy- 
bridization functions, but instead is dictated by the char- 
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Fig. 1: Schematic low-T phase diagram of a pseudogap Kondo 
impurity as function of chemical potential /i and Kondo cou- 
pling J. LM denotes the local-moment phase of an unscreened 
impurity, while SC and ASC are strong-coupling phases with 
screening. In the shaded regime (bounded by lines \p\ oc 
| J — J c \" , with v = 1 for r — 1) close to the critical coupling 
J c , the universal prediction Tk = «±|mI holds for r < 1. At the 
upper critical dimension r = 1, this is replaced by Tk = k\p\ 
(Tk oc \p\ x ) for the two signs of /i. Strong electron-hole asym- 
metry of Tk (/i) is present for a much larger range of parameters, 
for details see text. 



acter of the relevant critical fixed point (dubbed ACR) 
in the pseudogap Kondo model, which has been shown 
to have maximal p-h asymmetry near r = 1 [11]. Scale 
invariance of the critical fixed point, realized for r < 1, 
implies that the Kondo temperature Tk is proportional 
to the gate voltage, Tk = K±|mI; with a universal (and 
p-h asymmetric) prefactor k [14, 15], in a regime close to 
criticality shown in Fig. 1. Remarkably, we find that log- 
arithmic corrections at the upper critical dimension r = 1 
conspire such that Tk = n\p\ still holds for one sign of p 7 
while for the other sign Tk oc \p\ x with a universal expo- 
nent x w 2.6. Even more importantly, this strong asym- 
metry in T K between electron and hole doping extends 
over a large parameter range also away from criticality 
and should thus be easily detectable in graphene experi- 
ments. We also show that impurity tunneling spectra are 
sensitive to the sign of p as well. 

While a number of theory works on Kondo impu- 
rities and STM in graphene have appeared previously 
[12, 16-20], some of which also discussing aspects of gate- 
tuned Kondo physics [12, 16], the dependence of T K on the 
gate voltage has not been studied systematically, and the 
connection to quantum criticality has not been made pre- 
cise. Note here that the often employed methods of poor- 
man's scaling and slave bosons provide a qualitatively in- 
correct description of the QPT near r = 1 . 

Model. — Consider the standard spin-1/2 Kondo 
model, H = Hb + T^mp, with 

n b = X> c L c fc-> n imp = j^-^c CT (o) t ^ c ^(o) (i) 

ka aa' 

with a bath DOS p(ui) at site where c CT (0) = 
i. c ka- futji is the vector of Pauli matrices. Com- 
plications arising from orbital degeneracies will be dis- 



cussed towards the end of the paper. For the undoped 
case, we assume p(ui) = Po\io\ r for |u;| much smaller than 
some crossover scale Ap, while we make no assumptions 
about p for large energies below the band cutoff A (where 
p{oS) can be p-h asymmetric). A finite chemical poten- 
tial p shifts the minimum of p away from the Fermi level, 
p(uS) = po| w + p\ r i such that p > corresponds to elec- 
tron doping; we focus on \p\ <C Ap. For graphene r = 1, 
Ap ss 0.5 eV [12] and A ss 8cV. The low-energy p-h asym- 
metry of the impurity problem arises from asymmetries of 
both the bath and the impurity (the latter coming from 
additional potential scattering). 

Chemical-potential tuned crossover. For p ^ 0, 

the impurity spin will be screened below a temperature 
Tk(p) for any J. While J > J c results in a finite Tk as 
p — > 0, J <C J c yields an exponentially small Tk, lnTk oc 
— l/|/i| r - this follows from a weak-coupling RG treatment 
(i.e. poor man's scaling) . Naively, one would expect only a 
weak p-h asymmetry in T K (p) (because of the low-energy 
p-h symmetry of p(u>) at p = 0) and, in particular, T K (/i) 
to be minimal at p = for J > J c . Below we show that 
these expectations are incorrect. 

RG analysis. — We start with the theoretically most 
interesting case J = J c . Scale invariance of the critical sys- 
tem at p = entails T K = k\p\ where n(r) is a universal 
prefactor [15], because p scales as an energy and is a rele- 
vant perturbation at the critical fixed point. Wc shall show 
that this scaling prediction is indeed obeyed for r < 1. 
Notably, the prefactor is p-h asymmetric, T K = k±|m| for 
p ^ 0, for r* < r < 1 where the critical behavior is con- 
trolled by the ACR fixed point [7, 11]. 

The universal critical theory for the QPT near r = 1, 
which is not accessible from weak Kondo coupling, has 
been worked out in Refs. [10, 11]. It is the theory of a 
crossing of singlet and doublet impurity levels minimally 
coupled to conduction electrons, or, equivalently, a U = oo 
Anderson impurity model. Using the notation of Rcf. [11], 
its impurity part can be written as 

Himp = e \a)(a\ + V [|a)(e|c CT (0) + h.c] (2) 

where | f), | J,), and |e) represent the three allowed im- 
purity states. e is the tuning parameter ("mass") of the 
QPT, i.e., the (bare) energy difference between doublet 
and singlet states. The QPT occurs at some £rj = e c , with 
screening present for £o > e c . The hybridization Vq is 
marginal (relevant) at r = 1 (r < 1), i.e., and r = 1 plays 
the role of an upper critical dimension. Due to U = oo, 
the critical fixed point is maximally p-h asymmetric. Note 
that the pseudogap Anderson and Kondo models at any 
fixed r share the same universality class. 

Ref. [11] presented a field-theoretic RG analysis of the 
model (2) for r < 1, in a double expansion in Vq and 
(1 — r). Here we generalize the RG to include a finite p, re- 
stricting ourselves to the simpler momentum-shell scheme 
with ultraviolet (UV) cutoff D. We introduce dimension- 
less couplings v, e according to Vq = D^ r ^ 2 p 1 ^ 2 v and 
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Fig. 2: Schematic RG flow [11] for the U = oo Anderson model, 
(2), for jj, = and r* < r < 1. e and « are the impurity level 
energy and hybridization, respectively. The full (open) circles 
are stable (unstable) fixed points; the thick line separates the 
flows to the LM and ASC phases, for details see text. For 
r — > 1, the ACR fixed points moves towards e, v = 0. 



Eq = De for the hybridization and mass, respectively. The 
one-loop RG equations read 



dv 



1 



d\nD 

de 
d\nD 



r -v+ V -^F 1 { i i/D), 



2 2 

-e + v 2 eF 1 {fi/D)+v 2 F 2 {n/D) (3) 



with F ia {y) = |1 + y\ r ± 2|1 - y\ r . The last term in 
de/d In D simply describes the level shift due to the real 
part of the bath Green's function. For fi = 0, the RG 
equations (3) reduce to those of Rcf. [11], with the flow 
depicted in Fig. 2. They yield a critical fixed point at 
v* 2 = (1 - r)/3 and e* = -(1 - r)/(3r) for r < 1 - 
this is the ACR fixed point. The fixed points e = ±oo 
correspond to the screened and unscreened phases, ASC 
and LM, respectively. 

For small \i ^ 0, the flow at Eq = e c deviates from the 
critical fixed point on the scale \n\ [21]. For negative fi, 
the flow is directly into the screened (ASC) phase, e —> oo, 
implying Tk ~ i.e., k_ = 0(1). For a more detailed 
analysis, we can deduce Tk from the flow of m = e/v 2 . 
From Eqs. (3) one finds 



dm 
d\nD 



(4) 



with the critical-point value m* = —1/r. For fi =/= 0, 
we apply a two-stage RG procedure [21]: first, the cutoff 
D is reduced down to fx, following the flow equations for 
fx = 0. The second stage starts at D = fj, (and m = m* 
at J = J c ), and follows Eqs. (3,4). For simplicity, we may 
approximate here ^(y) ~ — \y\ r , which allows to integrate 
Eq. (4) analytically. Finally, Tk may be defined as the 
value of D where m has changed by unity, to(Tk) = m* + l, 
which yields 



T K = \p\[W{{r-l}/e)/{r-l)]^ 



(5) 



where W(x) is the Lambert W (or product log) function 
and e = 2.71828 . . . 

In contrast, for small positive n the flow at e c is driven 
towards negative e, i.e., the doublet (LM) fixed point, 



which implies that a spin 1/2 degree of freedom is formed 
on the scale This moment is subsequently screened 
due to the finite DOS at the Fermi level - this second 
crossover is akin to conventional Kondo screening. It in- 
volves restoration of p-h symmetry and is therefore not 
described by the RG equations (3). To obtain a rough es- 
timate of Tk here, we may again envision a second stage 
of RG, which starts at D = /i and now uses the language 
of poor man's scaling. This results in a conventional expo- 
nential expression for Tk, here to be used with bandwidth 
[i and dimensionless hybridization ~ \J\ — r (taken from 
the fixed-point value v*). Consequently, Tk = with 
ln«+ oc -1/(1 - r), i.e., T K < for r < 1. 

For r = 1, the upper critical dimension of the prob- 
lem, one may expect logarithmic corrections due to the 
logarithmic critical flow towards e = v = 0, e.g., v 2 = 
F 2 /[l - 3V 2 ln(T>/T> )] where D is the initial UV cut- 
off [11]. Remarkably, m* = —1/r also applies to this 
logarithmic flow, and hence logarithmic corrections are 
absent from the Tk expression (5) for fi < [22]. For 
fi > 0, the above two-stage RG analysis now suggests 
lnre + oc ]n(fx/Do) because the relevant dimensionless hy- 
bridization is v 2 cx — 1/ ln(/u/Do). Hence, n+ becomes 
scale-dependent: the scaling prediction is spoiled and re- 
placed by Tk/Dq cx \/j,/Dq\ x where x > 1 is now a univer- 
sal exponent. (To obtain an estimate for x here, a detailed 
analysis of the crossover between the two RG stages would 
be required which we do not attempt.) Tk oc \fj,\ x also im- 
plies that the low-energy impurity physics for fi > is 
characterized by two distinct scales, \x and Tk- 

Let us summarize the main physics at this point. The 
fundamental asymmetry between electron and hole doping 
arises from the p-h asymmetric character of the critical 
theory (2). Its phase transition can be driven by varying 
the real part of the bath Green's function. Now, varying 
/x leads to a linear variation of this real part, and hence 
drives the critical system either towards the singlet or the 
doublet phase on the scale We naturally expect that 
the electron-hole asymmetry, existing at the critical cou- 
pling, pertains to the off-critical situation as well. Our 
numerical results (Fig. 5) show that this is indeed the case. 

For completeness, we also mention the result for Tk(/u) 
at criticality for small r. Here, the slave-boson mean-field 
treatment of (1) [5] gives qualitatively correct results for 
the p-h symmetric critical (SCR) fixed point [7,11]. The 
mean-field equation determining Tk can be written as 



1 f , p(uj) , u> 

- = / efcji--^tanh 

J J w 2T K 



(6) 



At fi = 0, setting Tk = gives the critical coupling J c . 
For a system with p{ui) = po|^ ± p\ r and J = J c , Tk(/m) 
follows from the dimensionless equation 



dx 



x±l 



■ tanh ■ 



2T K 







(7) 



where x = to/\fi\, Tk = Tk/|/^|. For r < 1, the ultraviolet 
cutoff of the integral can be sent to infinity. Note that the 
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Fig. 4: Universal prefactors n of Tk = K±|/*|, obtained from 
NRG [24], for the symmetric (SCR) and asymmetric (ACR) 
critical fixed points. Both panels show the same data, with log 
and linear k axes, respectively. Tk is defined here via thermo- 
dynamic observables [26]. The lines represent the asymptotic 
forms exp( — 1/r) (solid), 1.8 exp[— 1.7/(1 — r)] (dashed), and 
3[W({r - l}/e)/(r - l)] 1/r (Eq. 5, dash-dot), for details see 
text. 



Fig. 3: NRG results [24] for impurity moment and entropy, 
Tximp and Simp, for a U = oo Anderson model (2) at the 
critical point, e c = —0.3105935, for different p, with dashed 
(solid) curves for p > (p < 0). The bath DOS is p(u) = 
\u+p,\Q(l— i.e., r — 1, and the hybridization V 2 = l/Vr. 
The screening process displays an extreme p-h asymmetry. 

non- universal factor po has dropped out. The result for 
Tk is independent of the sign of p, Tk = n\p\. Solving 
Eq. (7) yields n oc exp(— 1/r) for small r. The same result 
can be obtained using a two-stage version of poor man's 
scaling. In the first stage, p = 0, the RG equation for the 
dimcnsionlcss Kondo coupling j is dj/dhiD = —rj + j 2 
[11]. In the second RG stage, starting at the scale p, the 
fixed-point coupling j* = r simply enters the conventional 
exponential expression for Tk, i.e., Tk = \p\ exp(— 1/r). 

NRG results. — All the above considerations are well 
borne out by simulations using Wilson's numerical RG 
(NRG) technique [23]. In Fig. 3 we display NRG results for 
the impurity magnetic moment T^imp and entropy Si mp of 
a U = oo Anderson model with symmetric r = 1 DOS at 
the critical coupling, e = e c , for different p. At the critical 
fixed point T^imp = 1/6 and Simp = In 3 [10]. Various 
features are apparent in Fig. 3: The observables deviate 
from their critical value for T < \p\. More importantly, 
there is a striking asymmetry between the two signs of 
p, e.g., the magnetic moment T\ increases (decreases) for 
p > (p < 0) once T drops below \p\. This reflects 
the fact that the system is driven towards the doublet or 
singlet phase depending on the sign of p. Further, for 
p > a clear two-stage crossover is seen, where a spin- 1/2 
moment is formed below T ~ \p\, with T%i mp ~ 1/4 and 
S'imp ~ In 2, which is only screened at very low T. The 
power-law behavior, Tk(p) oc p x for p > 0, is visible, and 
a fit gives x ~ 2.6 ±0.1. 

Fig. 4 shows the universal prefactors n{r), obtained from 



NRG, of the "critical" relation Tk = ft±M- Note that 
Tk has been extracted from Tximp [26]; in general, Tk is 
only defined up to a prefactor of order unity. We have 
verified universality of k by performing calculations for 
both Kondo and U = oo Anderson models and for vari- 
ous power-law DOS with different high-energy parts. The 
results for n nicely follow the asymptotic forms derived 
above, i.e., appear well fitted by «(r — > 0) ~ exp(— 1/r), 
n + (r — > 1) oc exp[— 1.7/(1 — r)], and K_(r — > 1) oc 
[W(l + r)/(l+r-)] 1 /<\ NRG results away from the critical 
coupling will be discussed in connection with graphene in 
the next section. 

We finally note that the role of the two signs of p is 
interchanged if the sign of the model's p-h asymmetry is 
switched. Our quoted signs are valid for a U = oo Ander- 
son model, and also apply for a Co impurity on graphene, 
assuming a Kondo model and an orbital El configuration 
of Co [12], see below. 

Application to graphene. — For Co in the cen- 
ter of a graphene's hexagon, the hybridization function 
is strongly asymmetric on the scale of a few eV [12]. We 
have used a p(ui) which resembles [27] the tight-binding 
fit to the LDA DOS in Fig. 2 of Rcf. [12] (for Co in El 
configuration) to calculate Tk in a standard Kondo model 
as function of Kondo coupling J and chemical potential 
p. Results are displayed in Fig. 5. The strong electron- 
hole asymmetry is apparent, which continues to exist away 
from the critical coupling. At J = J c the linear and power- 
law behaviors of Tk(p) are nicely visible. Moreover, the 
minimum of Tk(p) for J > J c is found at significant posi- 
tive bias. 

A few comments are in order: (i) For simplicity, we have 
assumed a simple SU(2) spin 1/2 impurity. According to 
Ref. [12], spin-orbit coupling reduces the symmetry of the 
Co impurity problem from SU(4) to SU(2) on a scale of 
60 meV. This implies that our results are experimentally 
relevant only below a temperature of ~ 600 K [29] . More- 
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over, the Kondo coupling in our SU(2) model has to be 
taken essentially by a factor of 2 larger compared to the 
bare Kondo coupling of the Co impurity to obtain the 
same Tk (because the impurity degeneracy multiplies the 
Kondo coupling within the poor man's scaling equation.) 
Hence, we predict a J c of 2.2 eV for Co which approxi- 
mately matches [28] the estimate in Ref. [12]. (ii) The 
existing STM data, with T K ~ 15 K at \fx\ = 0.2 eV [4], 
can cither be consistent with our results for fi < and J 
significantly smaller than J c or with \x > and J « J c . 
The two cases differ in the change of Tk upon inverting 
the sign of li: in the former case Tk should become van- 
ishingly small, but huge in the latter (which would still 
be measurable if /i is reduced), (iii) The results in Fig. 5 
are switched according to/io —\i if the Co impurity is in 
an E2 configuration, because the hybridization functions 
for El and E2 symmetry are approximately identical after 
p-h transformation [12]. 

Tunneling spectra. — The different screening pro- 
cesses for fi ^ result in distinct spectral functions as well. 
In Fig. 6 we show NRG results for the impurity T matrix 
at T = 0, for situations with similar Tk, but reached with 
different combinations of fi and J. Most striking is again 
the case J = J c , where the high-energy part of the spec- 
trum follows the critical law ImT(w) oc 1/(cj| lnw| 2 ) [6,11]. 
For /j, < 0, the RG flow from critical to singlet behavior 
results in a large peak in ImT(w) of width Tk away from 
the Fermi level, located at u> ~ Tk ~ — /i. In contrast, the 
two-stage screening for fi > yields a more conventional 
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Fig. 5: NRG results for the Kondo temperature as function of 
fj, for different values of the Kondo coupling J, calculated for 
a DOS appropriate for Co on graphene [27] where J c w 4.3 eV 
[28]. Both top and bottom panels show the same data, with 
log and linear axes, respectively. 
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Fig. 6: T matrix spectral function for a Kondo impurity cou- 
pled to a DOS describing Co on graphene [27], for different 
values of the chemical potential and the Kondo coupling. All 
cases have a Tk between 20 and 30 K, but very different spec- 
tral shapes. The inset shows the same data on a larger scale. 
Note that the /i > spectra have a broad peak around (—(J,) 
in addition to the Kondo peak. 



Kondo peak at the Fermi level, but in addition a broad 
peak at a much larger energy w ~ —li (see inset of Fig. 6). 
This behavior persists away from J = J c , i.e., the Kondo 
peak has a pronounced asymmetry w.r.t. the Fermi level 
for li < (and J not too far from J c ). 

The impurity T matrix spectrum, Fig. 6, is in general 
not identical to the signal in a STM experiment, due to 
additional tunneling contributions into the host orbitals. 
The interference between the two is known to convert 
Kondo peaks into Fano lincshapes, with the Fano param- 
eter depending upon details of the tunneling process [30] . 
Nevertheless, the position and width of the Kondo feature 
are independent of this interference. Hence, we predict 
that for J ~ J c the Kondo feature in STM will either be 
centered at the Fermi level or be located at an energy of 
order Tk away from the Fermi level with a sign opposite 
to ll, depending on the sign of ll. 

Conclusions. — We have discussed Kondo screening 
for magnetic impurities in doped pscudogap Fermi sys- 
tems, with focus on the effects of nearby quantum crit- 
icality. For bath exponents r < 1, we found that the 
Kondo temperature follows the universal law Tk = 
for ll ^ in the critical regime. For r = 1, a subtle inter- 
play of p-h asymmetry and logarithmic corrections at the 
upper critical dimension leads to an extreme p-h asym- 
metry in Tk(/x), with linear or power-law behavior being 
realized depending on the sign of \x. 

For a situation corresponding to Co in graphene, we 
have predicted both the dependence of Tk on [i and char- 
acteristics of low-temperature STM spectra. Most striking 
is the strong asymmetry between electron and hole dop- 
ing over a large parameter range. This behavior is sur- 
prising because the relevant graphene DOS is weakly p-h 
asymmetric at low energies and hence a strong asymmetry 
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is not expected from standard weak-coupling approaches. 
Instead, it is the strong p-h asymmetry of the critical fixed 
point of the pscudogap Kondo model which is responsible 
for the physics advertised here (and which is not captured 
by slave-boson or weak-coupling approaches to the Kondo 
problem). 

Our findings call for systematic measurements of Xk as 
function of fi, in particular for gate voltages in the linear 
DOS regime, i.e., below 0.2eV. The asymmetry of Tk(h), 
together with the ab-initio data of Rcf. [12], will allow to 
determine the orbital nature of the Co impurity. More- 
over, a linear or power-law behavior of Tk(/i) will be a 
strong indication of nearby quantum criticality at /i = 0. 
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